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' stanton zero modes and study their symmetry properties. We present simula- 

O tions of the instanton ensemble for different numbers of Majorana fermions in 
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pose we calculate matrix elements of the adjoint Dirac operator between in- 



the adjoint representation. These simulations provide evidence that instan- 
tons lead to gluino condensation in supersymmetric gluodynamics. 

I. INTRODUCTION 

In order to improve our understanding of non-perturbative phenomena in QCD it is 
useful to view QCD from a larger perspective, as a member of a family of QCD-like theories 
with different matter contents. In this context we would like to understand the phase 
structure of QCD-like theories with Nf fermions in the fundamental representation of the 
gauge group and N a d fermions in the adjoint representation. Theories with adjoint fermions 
are special because the action may display a symmetry that connects bosonic and fermionic 
degrees of freedom, supersymmetry. Supersymmetry imposes powerful restrictions on the 
structure of the low energy effective action. These constraints have been used, for example, 
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to determine the phase structure of N = 1 supersymmetric QCD with N c colors and Nf 
flavors of fundamental quarks 

In addition to that, supersymmetry provides the opportunity to isolate certain non- 
perturbative effects, in particular instantons. This idea has been used in order to calculate 
the gluino condensate in the simplest supersymmetric gauge theory, SUSY gluodynamics. 
The strategy behind the so called weak coupling instanton (WCI) calculation f3)-[7| is to add 
to the theory a fundamental fermion together with its scalar superpartners and consider 
the regime where the expectation value of the scalar field is large. In this case there is a 
unique non-perturbative superpotential induced by instantons. Since the scalar vev is large, 
instantons are semi-classical and the superpotential can be calculated reliably. The superpo- 
tential determines the gluino condensate in the theory with additional matter. Finally, the 
matter fields can be decoupled by sending their mass to infinity. The result for the gluino 
condensate in SU{2) supersymmetric gluodynamics is 

(AA) = 6A 3 , (1) 

where A is the scale parameter defined in Q 

, „ /16ir 2 \ 1/3 / 8tc 2 \ 

h=Mpv [~w) exp rw (2) 

Here, b = 3N C is the first coefficient of the beta function in supersymmetric gluodynamics 
and Mpv is a Pauli-Vilars regulator. 

There is an old puzzle concerning this result. The puzzle is connected with the fact that 
there is an alternative method for calculating the gluino condensate, usually referred to as 
the strong coupling instanton (SCI) method |@-12|, for a review see ||13|| . In the case of N c = 2 



one considers a four fermion correlation function. This correlation function is a topological 
quantity. Not only can it be saturated with one instanton, but supersymmetry implies 
that the correlator is just a constant. At short distance, one expects that this constant is 
saturated by small instantons and can be calculated reliably. The gluino condensate is then 
extracted by using clustering. The puzzle is that the result differs from the WCI calculation 
by a factor 4/5. 



Several suggestions have been put forward in order to resolve the puzzle []T4| — pTTfl . We 
do not wish to discuss these possibilities in detail. Instead, we would like to employ a 
somewhat different, more qualitative approach. Even though there is no direct instanton 
contribution to the gluino condensate, one would still expect configurations with instantons 
and anti-instantons to contribute to the gluino condensate. Here we have in mind that 
the theory is studied in a finite volume and in the presence of a non-zero mass term. The 
thermodynamic limit is approached by taking the volume to infinity before we let the mass 
go to zero. The mechanism for gluino condensation is similar to the instanton liquid picture 
of quark condensation in ordinary QCD [|18| -|20|1. For Nf > 1 there is no direct instanton 



contribution to the quark condensate but chiral symmetry breaking may take place in an 
ensemble of instantons and anti-instantons in the thermodynamic limit. The Banks-Casher 
relation (qq) = —np(0) |^TJ connects the quark condensate with the density of eigenvalues 
of the Dirac operator at zero virtuality. For simplicity let us consider an ensemble with an 
equal number of instantons and anti-instantons. In this case the Dirac operator no longer 
has any exact zero modes. However, if the interaction between the instantons is sufficiently 
weak, the approximate zero modes associated with individual instantons and anti-instantons 
form a zone around zero virtuality and lead to spontaneous chiral symmetry breaking. This 
quark condensation mechanism has been investigated in some detail, both analytically and 
on the lattice |2^j23|1 , and the results seem to support the instanton picture. 



In the present work we wish to extend these studies to theories with fermions in the 
adjoint representation. Since we are dealing with a strongly coupled theory, our calculations 
are necessarily approximate. In particular, we will have to restrict ourselves to the contri- 
bution of small instantons for which the semi-classical description is appropriate. On the 
other hand, the methods we are using are applicable also to non-supersymmetric theories 
with several flavors of adjoint fermions. In addition to that, we can use these methods to 
study non-constant correlation functions that determine the spectrum of the theory. 

The paper is organized as follows. In section |TJ we discuss some general aspects of chiral 



symmetry breaking in theories with fermions in the adjoint representation. In section fl I 



we describe the structure of the instanton zero mode wave functions and in |V| we calculate 
matrix elements of the Dirac operator between zero mode states. These results are used 
in order to determine the fermion determinant in the field of an instanton-anti-instanton 
pair (section [V]) and to calculate the gluino condensate in a random instanton ensemble 
(section |V|). In section |V11| we describe simulations of an interacting instanton ensemble 



with different numbers of fermions in the fundamental and adjoint representation. 



II. QCD WITH ADJOINT FERMIONS 

QCD with adjoint fermions is defined by the lagrangian 

N ad 1 _ 1 

£ = E ^2W 6 Ag 6 - —G%G%, (3) 

i=i z 4 y 

where X a M is a Majorana fermion in the adjoint representation of the gauge group and G* 
is the usual field strength tensor. The covariant derivative in the adjoint representation is 
given by 

Df = d^5 ab + f abc A^ (4) 

For several Majorana flavors the theory (§) possesses a SU{N cu j) chiral symmetry. A non-zero 
gluino condensate 

(A2A$> = -Pa (5) 

breaks this symmetry to SO(N a d) p4j. This fact can be seen most easily by considering the 
conserved vector and axial- vector currents p5| . There are \N a d{N a d — 1) conserved vector 
currents V' l J = A^T/jA^ and \N a d{N a d + l) classically conserved axial- vector currents A l j = 
Am7^75Am- The singlet axial current A" is anomalous. At the quantum level this leaves 
N% d — 1 conserved charges that generate the SU(N a d) chiral symmetry. Gluino condensation 
breaks the axial symmetries and leads to the appearance of \{N"^ d + N a d — 2) Goldstone 
bosons. The unbroken \N a d(N a d — 1) vector charges generate the residual 0(N a d) symmetry. 
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In the case of supersymmetric gluodynamics, N a d = 1, there is no continuous symmetry. 
Instantons break the axial U(1)a symmetry but leave a discrete Z^ c symmetry intact. This 
discrete symmetry is spontaneously broken by gluino condensation. As discussed above, 
the value of the gluino condensate is known from a weak coupling instanton calculation. 
There are no predictions for the spectrum of the theory, but we expect the lowest states to 
fill out a chiral supermultiplet containing a scalar and a pseudoscalar meson as well as a 



Majorana fermion. These results can be summarized in terms of an effective lagrangian p6 
This is not an effective lagrangian in the Wilsonean sense. The effective action does not 
generate the low momentum scattering amplitudes of the theory. Instead, it mainly serves 
as a generating functional for the anomalous Ward identities of the theory. 

III. INSTANTON GAUGE POTENTIAL AND FERMIONIC ZERO MODES 



In theories with adjoint fermions it is convenient to employ a spinor notation for spin, 
vector, and color indices ||. We can convert vectors to spinors using 

V a a = V^a^aa- (6) 

The euclidean spinor conventions used in this paper are summarized in Appendix A. The 
instanton gauge potential couples spin to color degrees of freedom. A field A a in the adjoint 
representation of 577(2) can be represented by a symmetric tensor A a @ 

A a = A^e^T*)}. (7) 

In spinor notation, the instanton gauge potential in regular gauge is given by 

A$ = -2i{S*aj + *M)-f^r (8) 
We can transform the gauge potential to singular gauge using the gauge transformation 

U« a = x,(a^f a . (9) 
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Note that this matrix transforms an undotted color index into a dotted one. We can perform 
a 'fake' conversion of the dotted spinor back to an undotted one using the fact that (a ) aa 
is just the unit matrix. 

In the case of SU(2), the Dirac operator in the background field of an instanton has four 
zero modes. The first two are conventionally called the supersymmetric (ss) zero modes [T3| 

a^ ) = (^ + «)7(^47F' (10) 

where (3 = 1,2 enumerates the zero modes. The other two are referred to as the supercon- 
formal (sc) zero modes 

A 75 . = (Px 6 - + 5 s x' 1 -)— - fill 

In singular gauge, the zero modes are given by 

A 2») = ^ + ^7S(?T^ < sc >' < 13) 

Analogously, we can construct the zero modes of the Dirac operator in the background 
field of an anti-instanton. The regular gauge supersymmetric zero mode has the structure 
X f) ~(^f + a*$,etc. 

The effect of the zero modes on the propagation of fermions can be summarized in terms 
of an effective lagrangian ||27J. The 't Hooft effective interaction in the case of one Majorana 
fermion in the adjoint representation of SU(2) was determined in ERKDI]. The result is 



C = ^ (-) 4 exp (--) p'dp {\ a M \ a M d^» M d»\» M + X a M l 5 KfdAl^X 



b 

M 



-AWWV^}. (14) 

This result has to be interpreted with some care. The notion of an effective interaction 
induced by instantons of some fixed size is incompatible with supersymmetry. In order to 
derive manifestly supersymmetric results we always have to integrate over the collective co- 
ordinates of the instanton. Nevertheless, it is instructive to compare the result (|i~4]) with the 
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effective interaction in the case of Nf = 2 Dirac fermions in the fundamental representation. 
The structure of the two interactions is quite similar, which suggests that instantons may 
lead to similar physical effects. The most important difference between the two effective 
lagrangians is the presence of derivatives acting on two of the four Majorana spinors in 
This difference is connected with the asymptotic behavior of the supersymmetric zero 
modes, which is not ~ 1/z 3 , but l/z 4 . 



IV. MATRIX ELEMENTS OF THE DIRAC OPERATOR 



In the following, we wish to study the spectrum of the Dirac operator in an instanton 
ensemble. For this purpose, we have to calculate matrix elements of the Dirac operator 
between the zero modes of individual instantons and anti-instantons 



t ia = i d A x\wm ab ^ b A- 



(15) 



An ensemble of instantons and anti-instantons is only an approximate saddle point of the 
action. If the system is sufficiently dilute then the instantons and anti-instantons are well 
separated and the approximate saddle point solution for the gauge potential is given by a 
simple sum of the gauge potentials of the individual instantons. For this purpose, the gauge 
potential of the individual instantons has to be put in singular gauge. In the sum ansatz, 
we can use the equations of motion of the fermion fields in order to replace the covariant 



derivative in (|15|) by an ordinary derivative 

T, , = - / <l''xX;(iil))X\ 



(16) 



The structure of the Dirac operator is dictated by the form of the zero modes. In the 
background field of an instanton-anti-instanton pair we have 



I A 







rpSS — SS rpsS — SC 

1 IA 1 IA 



rpSC—SS rpSC—SC 

1 IA 1 IA 



\ 



T 



Al 



rpsS — SC 

1 AI 



T 



SC—SS rflSC—SC 



Al 



T 



Al 



(17) 



where the matrix elements T|j, . . . are real quaternions. These quaternions can be decom- 
posed as 

(nrv = T "KV/3 ( 18 ) 

(rpsS — SC\ rpsS—SC i rpsS — SC (— \7 ^Ofl^ 

l J AT J/3/3' — J e /3/3' "T" J a*i/ e p-y\ <J p.v) pn l zu J 

(?ir s W = T sc ~ ss eppt + r;r ss (^)/6 7 ^, (21) 

Here, T^ s and T^ c are real vectors, T ss ~ sc and T sc ~ ss are real scalars and T^~ sc and T^~ ss are 
self-dual and anti-self-dual tensors, respectively. Chiral symmetry implies that the diagonal 
blocks of are zero. The upper right and lower left blocks are related by hermitean 
conjugation. For example, we find that 

(7>T SS W = T?{<r/P. (22) 

In general, we have (T')ai = (T)ia- The eigenvalues of fllTf) come in quartets (£,£, — £, — £). 
These results are in agreement with the general arguments presented in PS| , PD[ . 

The functions T^ s , . . . depend on the collective coordinates of the instanton and anti- 
instanton. We will characterize the relative color orientation by the four vector w M = 1/2 - 
tr(U A Uja^). Here, Uj t A are SU(N C ) matrices that characterize the color orientation of the 
instanton and anti- instanton. For color SU(2) is a real vector with u 2 — 1. Using 
rotational symmetry and the fact that Tai is quadratic in we have 

=vTf + u (1 (u-z)TZ a + z,(u-z) 2 T°\ (23) 
T° c =z,Tr + u,(u-z)T° c + z ll (u-z) 2 n c , (24) 

TSS _ SC = T ss-sc + T ss-sc( u . ^ ( 25 ) 

sc = («^ - uj,) (u ■ z)Tr sc , m) 

where z il = — z 1 ^ and the functions T( s , . . . depend on (|^|, pi, Pa)- For simplicity, we will 
assume that the dependence on pj tA only enters through their geometric mean p = ^fpipA- 



The fact that this assumption is valid to fairly good accuracy was checked in the case of 
fundamental fermions. In a more sophisticated treatment of the instanton-anti-instanton 
gauge configuration the dependence of the overlap matrix element on (z, pi, pa) is restricted 
by conformal invariance ||31|| . 

In the following we shall outline the calculation of the invariant functions T£ s , .... We 
describe the case Tf* s in some detail but relegate the results for the other functions to 
appendix B. Using the expression ( 12"1) for the wave function of the supersymmetric zero 
modes in singular gauge we find 

T° s = jtr (o^apapOo) tr {a v o a o^o^) + tr (p p a v a a o a ) j 

■u p u a I (fx(t)p V (x - z)d a (t)p 1 (x), (27) 



where <fip V {x) is the profile function of the supersymmetric zero mode 

M*) = ~ , ffW (28) 

7T [X 2 + p A Y 

The integral in ( P7|) is most easily calculated in Fourier space. The Fourier transform of tp^ 
is given by 

<f>tu,{k) = SpuMk) + kpk u <p 2 {k) (29) 

with 

2np 2 [4 / 4 \ r ^ 2 



( fc )= -f" {-,- ) K x {y) - -K G {y) , (30) 



y [y 3 \y 2 ) y 



Mk) = -2.pl - (i? + lj Ki(y) - (1 + ij K (y)} , (31) 

and y = fcp. K n (y) is the modified Bessel function of the first kind and order n. We can 
now calculate the overlap integral and perform the traces. In momentum space the result is 
given by 

T**(k) = H) (-2k v - 8u v (u ■ k) + I6(u ■ k) 2 ) \Mk)\ 2 - (32) 
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Finally, we can determine the functions T( s 2 3 by performing the inverse Fourier transform. 
In the d 4 k integral all integrations except for the one over the absolute magnitude of k can 
be performed analytically. We find 

= ^-2j dk \^k{hz) - 16^j 2 (*s) j \Mk)\ 2 (33) 
T ?( z ) = ^J dk l^hikz) - 32^j 2 (A*)j \Mk)\ 2 (34) 
T i S (z) = ^J dk Wk 4 h(kz) \Mk)\\ (35) 



where j n (x) is the spherical Bessel function of order n. The integrals (33) have to be 
performed numerically. The results are shown in Figure 1. In the following, we will use a 
simple parameterization of the numerical results. In the case of the supersymmetric overlaps, 
we use 

—1.265 

pT ' S ^ = 1.0 + 2.34z 2 + 0.35z 4 + 0.24*6' (36) 

- pT ss (z) L°5£ -6.36** (37) 

9 2 [ ] ~ (1.0 + 0.38^)3 + (i.o + 0.68z 2 ) 4 ' 1 ' 

*™ = (lo + oIL-t ' (38) 

where z = z/p. These parameterizations respect the asymptotic behavior of the overlap 
integrals. In particular, we have T ss (z) ~ 1/-2 5 , T sc (z) ~ 1/z 3 and T ss ~ sc (z) ~ l/z A . 

For completeness, let us compare these results to the corresponding expressions in the 
case of fundamental fermions. In this case, there is only one zero mode per instanton. 
The overlap matrix element Tja is a real number in the case of SU(2), and complex for 
SU(N C > 2). Tia satisfies the symmetry relation Tia = T* Al . As a consequence, the 
eigenvalues are real and occur in pairs (£, — £). We can extract the dependence of Tja on 
the collective coordinates. The result is 

1 Az 

T AT d = 0- z)T f (z, Pl ,p A ), T f (z, Pl ,p A )^ , 2/ , Y\2 ' ( 39 ) 

P/Pa (2.0 + ^/(p/Pa)) 2 
We note that the fundamental overlap matrix element only depends on one SU(2) angle 
cos# = (u ■ z). From the asymptotic form of the zero mode solution one finds T* (z) ~ 1/z 3 . 
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V. THE FERMION DETERMINANT IN THE FIELD OF AN 
INSTANTON-ANTI-INSTANTON PAIR 



Before we study gluino condensation we would like to make a brief digression and discuss 
the gluino induced instanton-anti-instanton interaction. This interaction will play an impor- 
tant role in the calculation of the gluino condensate in an interacting instanton ensemble. 

The probability to find an instanton-anti-instanton pair characterized by the collective 
coordinates (zi t A, Pi,a,Ui,a) is controlled by the weight factor exp(— S) det(i^ + m). The 
first factor is the well known gluonic interaction. If the instanton and anti-instanton are 
well separated it has the dipole form p2| 



S = 2S - (l-4 cos 2 6) , (40) 



2 4 



where S = (8n 2 )/g 2 is the single instanton action and cos# is the SU(2) angle introduced 
above. We note that the interaction is attractive if the color orientation is aligned with the 
spatial orientation, cos 6* = ±1. The second factor is the fermion determinant. In the case 
of fundamental fermions, it is also well known. We have 

16 z 2 

det(O) = cos 2 9^- j, (41) 

W> p 2 I p A { < 2.Q + z 2 /{p I p A )f ^ ' 

which is also attractive for cos# = ±1. We also note that the interaction peaks at z 2 ~ PiPa- 
Using the results of the last section we can calculate the fermion induced interaction 
with fermions in the adjoint representation. In order to calculate the determinant for one 
Majorana fermion we take the square root of the corresponding expression for a Dirac fermion 
in the adjoint representation. For simplicity, let us begin with the determinant in the basis 
of the supersymmetric zero modes only. We find 

det(^) ss = I (Tf) 2 + ((Tf) 2 + 2T?T? + 2T 1 SS T 3 SS ) cos 2 6 (42) 

+ ((T 3 SS ) 2 + 2TfT 3 ss ) cos 4 9 . 



The result for the superconformal zero modes is even more simple, 
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det(^) sc = ((Tf) 2 + ((Tf) 2 + 2Tl c T° c ) cos 2 o\ . (43) 

This expression is quite similar to the determinant for fundamental fermions. The super- 
symmetric determinant ( fH|) is somewhat more complicated, but also peaked for cos# = ±1. 
When the mixing between supersymmetric and superconformal zero modes is included the 
fermion determinant depends on other SU(2) angles in addition to cos 8. We show nu- 
merical results for log(det(i^)) as a function of z, cos# and cos0 in Fig.^. Here we have 
taken z M = zb^ and defined cos# = w 4 and sin 6 cos = u 2 - We observe that again the 
determinant peaks for z 2 ~ PiPa- For large z, the determinant behaves as ~ 1/z 16 . More 
importantly, we find that the interaction is again most attractive for cos# = ±1. There is 
some dependence on cos0, but it is not as pronounced as the dependence on cos 6. This 
means that the gluino induced interaction for one Majorana fermion is qualitatively similar 
to the quark induced interaction with an effective number of quark flavors between Nf = 2 
(which gives det ~ 1/-2 12 ) and Nf = 3 (corresponding to det ~ 1/z 18 ). 

VI. GLUINO CONDENSATION IN A RANDOM INSTANTON ENSEMBLE 



In this section we study gluino condensation in a random instanton ensemble. This means 
that we will assume that the collective coordinates of the instantons and anti-instanton are 
distributed randomly. In particular, we shall neglect the effect of the fermion determinant 
on the distribution of instantons. This is not a good approximation even in ordinary QCD 
and it certainly cannot be correct in a supersymmetric theory. Nevertheless, using the 
approximation of randomness we can get some analytic understanding of the dependence of 
the gluino condensate on the parameters characterizing the instanton liquid. We can also 
get an estimate of the relative size of the quark and gluino condensates in theories with both 
fundamental and adjoint fermions. 

The simplest model of the spectrum of the Dirac operator is based on the assumption 
that the non-zero matrix elements of the Dirac operator are Gaussian random numbers ||T9 
The distribution is characterized by the first moment 
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^ 2 =(|tr(TtT)). (44) 

The eigenvalue distribution for the Gaussian ensemble is given by a semi-circle where the 
density of eigenvalues at zero virtuality is p(0) = (N/V)(na)^ 1 . Here, (N/V) is the number 
of eigenstates per unit volume. The first moment of the overlap matrix can be estimated 
by averaging | Ta/ | 2 over the collective coordinates of the instantons. Using ( |3~9"1) we find the 
first moment of the Dirac operator for fermions in the fundamental representation of SU(2) 

m 



where p is the average size of the instanton. This parameter, just like the density of in- 
stantons, cannot be determined in the semi-classical approximation. In the instanton liquid 
model of the QCD vacuum it is assumed that p = 1/3 fm and (N/V) = lfm 4 [|18[]. Using 
these values we find 

i /N\ 1 / 2 

(qq) = --3 1 / 2 ( ) ~ -(230 MeV) 3 , (46) 



np V V . 

in very good agreement with the phenomenological value (which, of course, applies to color 
ST7(3)). 

The same arguments can be applied to gluino condensation in a random instanton en- 
semble. In this case we need to determine the first moment of a quaternionic matrix with 



the matrix elements determined in section IV. We find 



o aA = [y) 0.43pvr, (47) 

which is somewhat smaller than (|45|). There are four times as many eigenstates per unit 
volume but for a Majorana fermion the Banks-Casher relation has an additional factor 1/2, 
(AA) = —tt/2 ■ p(0). We finally get the following estimate of the gluino condensate 

1 /N\ 1/2 

(XX) = — —4.6 ( ) . (48) 



np \V 



13 



Here and in what follows we have dropped the subscript M on the Majorana spinors. This 
result is a little more than twice as large as the corresponding result for a Dirac fermion in 
the fundamental representation. As we saw, this is mainly due to the effective number of 
zero modes in both cases. We emphasize that the gluino condensate is proportional to the 
square root of the instanton density, which is also what one would expect if the condensate 
is extracted from the four-point function using clustering. 

We have checked these estimates by performing a numerical calculation of the spectrum 
of the Dirac operator in a random instanton ensemble. This means that instead of assuming 
the matrix elements of the Dirac operator to be random we take the collective coordinates 
of the instantons and anti-instantons to be random. We calculate the spectrum of the Dirac 
operator and determine the gluino condensate using 

< AA M/^ (A) OT (49) 

The results are shown in Fig. [| We observe that the spectrum is not a semi-circle but is 
peaked towards zero virtuality [Q. This non-analyticity is smoothed out when we calculate 
the condensate for a non-zero quark or gluino mass. Again using (N/V) = lfm 4 and 
p = 1/3 fm as well as m q = m g = 20 MeV we find (ipip) = -(260 MeV) 3 and (AA) = 
-(347 MeV) 3 . 

VII. GLUINO CONDENSATION IN AN UNQUENCHED INSTANTON 

ENSEMBLE 

As we stressed in the previous section, the assumption of randomness is not expected 
to be very useful. The fermion determinant is given by the product of all eigenvalues of 
the Dirac operator, while the quark or gluino condensate is determined by the density of 
small eigenvalues. This implies that the determinant tends to suppress fermion condensates. 
In particular, we expect that the strength of chiral symmetry breaking is reduced as the 
number of fermion flavors is increased. 
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In this section we shall study this problem using simulations of the instanton ensemble 
in QCD with fundamental and adjoint fermions. We consider the partition function 

Z = J (li d ^d( Pi )^ deb(ty f + m q ff det(# a + m 9 ) N ^ /2 exp(-S) (50) 

Here, Q denotes the collective coordinates of the instanton, d(p) is the single instanton 
distribution |27|j29|j3^1 , tff, a are the Dirac operators in the fundamental and adjoint repre- 



sentation, and exp(— S) is the gluonic interaction between instantons. In order to study 
spontaneous symmetry breaking in a finite volume we introduce non-zero quark and gluino 
masses m q , g - We will study the limit m q , g — > in some detail. 

The partition function ( |50| ) suffers from the usual IR problem connected with large 
instantons for which the semi-classical approximation does not apply. In practice, we deal 
with this problem by introducing a short range repulsive core in the gluonic instanton 



interaction, see section V.C. in |22| for a more detailed discussion. The repulsive core 
eliminates the contributions of large instantons and very close pairs. This particular method 
for suppressing objects that are not semi-classical has the virtue that it respects the classical 
scale invariance of Yang-Mills theory. 

The instanton ensemble is characterized by two numbers, the scale parameter A that en- 
ters into the instanton weight d(p) and a dimensionless parameter A which determines the 
size of the core. Lacking a better theory of topological fluctuations beyond the semi-classical 
domain we have to fix A phenomenologically. This could be done, for example, as soon as 
lattice information on the spectrum and other properties of theories with adjoint fermions 



becomes available [^,^]. In this work we will use the same value that was employed in 



studies of QCD with fundamental fermions. It leads to a dilute instanton ensemble charac- 
terized by the dimensionless parameter p 4 (N/V) ~ 0.12. For simplicity we will concentrate 
on simulations at a fixed instanton density (N/V) = 1.0 A 4 . 

To set the stage, we show results for Nf = 1, . . . , 4 flavors of fundamental fermions. Fig. 
[| shows the quark condensate as a function of the quark mass from simulations in a euclidean 
box of size V = 2.0 4 A 4 . The case of only one flavor is special. The chiral condensate persists 
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even if the limit m q — > is taken in a finite volume. This is due to the fact that for Nf = 1 
the quark condensate is dominated by direct instanton contributions. The result for Nf = 2 
is characteristic of spontaneous symmetry breaking. The quark condensate vanishes as the 
quark mass goes to zero but shows a clear plateau for larger quark masses. One can verify 
that the onset of chiral symmetry breaking moves towards smaller masses as the volume is 
increased. For more than two flavors the chiral condensate is significantly reduced. In the 
case of three flavors the signal is already quite weak. Using simulations in bigger volumes 
one can verify that chiral symmetry is indeed broken. There is no clear evidence for chiral 
symmetry breaking in simulations with four or more flavors. 

Fig. |]b shows the gluino condensate measured in simulations with one or two flavors 
of Majorana fermions in the adjoint representation. For N a ^ = 1 there is clear evidence 
for spontaneous symmetry breaking. Indeed, the behavior is more reminiscent of the case 
Nf = 1, where (qq) receives direct instanton contributions, than the case Nf = 2, in which 
chiral symmetry breaking is a collective effect. 

These observations can be understood in more physical terms. Supersymmetric gluody- 
namics has no Goldstone bosons, so finite volume effects are much weaker than in Nf = 2 
non-supersymmetric QCD. This means that in a fixed volume, gluino condensation can be 
observed for gluino masses that are significantly smaller than the quark masses required to 
produce quark condensation. In the standard picture, there is a discrete chiral symmetry 
which is broken by gluino condensation. This means that if the gluino mass is too small then 
chiral symmetry will be restored because of tunneling between the Z 2 vacua. This is dif- 
ferent from Nf = 1 non-supersymmetric QCD where instantons leave no unbroken discrete 
symmetries. 

The value of the gluino condensate is (AA) ~ 2A 3 . This result has the correct order of 
magnitude but it cannot yet be compared directly to the prediction ([!]). First of all, we 
use a different definition of the scale parameter. In order to make contact with our work 
on QCD we use a Pauli-Vilars scale parameter. Second, we have an additional parameter 
A which controls the boundary of the semi-classical regime. Finally, we have performed 
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the simulations at a fixed density of instantons (N/V) = 1.0 A 4 . It is this choice which 
effectively sets the scale in our calculation. 

In Fig. |]b we also show the gluino condensate measured in simulations with N a d = 2 
Majorana flavors. The condensate is very small and there is no clear evidence for spontaneous 
chiral symmetry breaking. 

The spectrum of the Dirac operator for Nf = 2 quark flavors and N a d = 1 Majorana flavor 



is shown in Fig. ||a and b. The spectra were determined in simulations with m qg = 0.1A 



-l 



Again, we observe that in both cases there is a finite density of eigenvalues as A — > 0. For 
Nf = 2 the spectral density near A = is flatfj, whereas in the case N a d = 1 it is growing 
towards small A. Again, this is similar to the case of only one fundamental fermion. The 
results are consistent with the effective field theory prediction [p?| , p8 



'(\-q) - S o (N f -2)(N f + P) 

Here, (3 is the Dyson index of the random matrix ensemble with the appropriate symmetry. 
We have (3 = 1 for fundamental fermions in SU(2), (3 = 2 for fundamental fermions in 
SU(N > 2), and (3 = 4 in the case of fermions in the adjoint representation. Nf denotes the 
number of Dirac or Majorana flavors in the cases (3 = 1,2 and (3 = 4, respectively. Sq is the 



magnitude of the quark condensate and /„. the pion decay constant. The expression (pi ) 
summarizes the fact that the spectrum is peaked towards small virtuality for both Nf = 1 
and N ac i — 1 while it is flat for Nf = 2. Effective field theory predicts the slope of the Dirac 
spectrum under the assumption that chiral symmetry is broken. The theory cannot predict 
whether chiral symmetry breaking takes place for some given Nf or N a d- 



1 Fig. H shows that the spectral density is flat for intermediate values of A. There is a finite 
volume suppression of the spectral density for small A and a 0(m 2 ) peak at A = 0. To show that 
the spectral density is flat at A = in the limit V — > oo, m — > requires more numerical work. 
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VIII. CONCLUSIONS 



In summary we have studied gluino condensation and the spectrum of the Dirac operator 
in an instanton ensemble. We employ the semi-classical approximation and focus on the 
Dirac operator in the subspace spanned by the zero modes of the individual instantons and 
anti-instantons. We have shown how the quaternionic structure of the Dirac operator in 
theories with adjoint fermions emerges naturally from the spin and color structure of the 
zero modes. The dependence of the matrix elements on the collective coordinates of the 
instantons is quite complicated but qualitatively similar to the simpler case of fundamental 
fermions. 

We have provided evidence that gluino condensation does take place in an ensemble of 
instanton and anti-instantons. In a random ensemble, the gluino condensate is proportional 
to the square root of the instanton density. In supersymmetric gluodynamics we find that 
gluino condensation persists even if interactions between the instantons are taken into ac- 
count. We observed that finite volume effects are much weaker than in QCD with two flavors 
of fundamental fermions. This is consistent with the fact that supersymmetric gluodynamics 
has a large mass gap. In QCD with more than one adjoint flavor we found no compelling 
evidence for gluino condensation. 

There are many problems that remain to be studied. In particular, it would be interesting 
to make a systematic study of gluino and gluino-glueball correlation functions. There are two 
types of correlation functions: Constant correlators that provide information on condensates, 
and x-dependent correlators related to the spectrum. These correlation functions will also 
show to what extent supersymmetry is realized in the limit m g — > 0. In addition to that, it 
would be interesting to search for evidence of Zi domains and to investigate the dependence 
of the results on the topological sector of the theory. In this work we have used the zero mode 
wave functions that correspond to trivial holonomy and anti-periodic boundary conditions 
on the fermions. This suggests the question of how the results are changed if the boundary 



conditions are modified. In this case, the zero modes discussed in [[!9[] will come into play. 
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Finally, it is important to study the role of very large instantons that were excluded in the 
present study 

There have been suggestions that objects with fractional topological charge may play a 
role in theories with adjoint fermions p!)|- ^ . |T5[ | . These objects can give a direct contribution 
to the gluino condensate. Because of tunneling between the different phases the presence 
of such objects cannot be inferred from the behavior of the gluino condensate as a function 
of the quark mass in a finite volume. One should be able, however, to detect the presence 
of fractionally charged objects in lattice simulations by looking for zero modes of the Dirac 
operator that do not appear in multiples of 2N C ||43|| . In this context it would also be 
interesting to study gluino condensation for N c > 2. For adjoint fermions the number of zero 
modes per topological charge increases with N c . One might therefore doubt that instantons 
alone are sufficient to trigger gluino condensation in large N c SUSY gluodynamics. It has also 
been suggested that fractionally charged objects can be thought of as instanton constituents 



^,|4T|f45||. One might then envision a situation where if N c is small, or instantons are 



small, fractionally charged objects are bound into instantons while for large N c , or for large 
instantons, topological objects dissociate and the instanton liquid should be replaced by 
liquid of fractional charges. 
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APPENDIX A: EUCLIDEAN SPINOR CONVENTIONS 



We use the following euclidean spinor conventions 



7m 



a, 



75 



'-10^ 



V 01 / 



a M ={ia,l), a li = (-ia,l), 



Mr i 



a/3 



Indices are raised and lowered with e a/3 and e alS where e^e^ = 8™ and e alS 
euclidean sigma matrices have the following properties 



2 £ pa® pa ■ 



(Al) 

(A2) 

(A3) 
(A4) 

e a/3 . The 

(A5) 
(A6) 
(A7) 
(A8) 



APPENDIX B: MATRIX ELEMENTS 

In this appendix we collect the remaining matrix elements of the Dirac operator. We 
define the profile function of the superconformal zero mode 



4>p 



X, 



V2n (x 2 + p 2 ) 2 
and its Fourier transform = —ik^^k) with 

Mk) = -V2np 2 K l (kp). 



(Bl) 



(B2) 



The matrix elements of the Dirac operator between superconformal zero modes are deter- 
mined by 
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T° c (k) = (+z) (-2k v + 8u v (u ■ kf) k\fa{k)\' 



(B3) 



and the matrix elements between supersymmetric and superconformal zero modes lead to 



T ss - sc {k) = (2 - 8(u ■ k) 2 ) kfaWMk), 

(k) = 8 (ujzy - uJkfA (u ■ k)k(j) 2 (k)(j)3(k). 



n-iss — sc 



(B4) 
(B5) 



From these results we can extract the invariant functions 



dk2k i j 1 {kz)\<j )3 (k)\' 



(B6) 



and T% c (z) 



-4T( c {z) as well as T 3 sc (,z) = 0. Also 



57T 

1 
^2 



A- 3 



2k A j (kz) - 8— j x {kz) 



dk 

dk8k 4 j 2 (kz)4> 2 (k)fo(k) 



(B7) 
(B8) 



and T| s sc {z) = — T 2 SS sc (z). Numerical results for these functions are shown in Fig. [I]. The 



results can be parametrized as 



-0.25,2 



1.0 + 0.422 2 + 0.2H 4 



(B9) 



as well as 



-0.17 



prr sc {z) 



1.0 + 0.05* 2 + 0.08z 4 ' 

l.2z 2 0.014^2 



(B10) 
(Bll) 



(1.0 + 0.45^2)3 (1.0 + 0.21z 2 )3 

where z = z/p. The overlap matrix elements T ss ~ sc are related to the corresponding 
functions with the supersymmetric and superconformal zero modes interchanged. We find 
T( c 2 ~ ss = -T? s f sc and T 3 S 



"ISC — ss 

3 



rrisS — SC 

J 3 
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FIGURES 

FIG. 1. Invariant functions characterizing the overlap matrix elements of the Dirac oper- 
ator. Fig. a,b and c show the diagonal overlap matrix elements between supersymmetric and 
superconformal zero modes, and the mixed supersymmetric-superconformal matrix elements. 

FIG. 2. Logarithm of the adjoint fermion determinant in the field of an instan- 
ton-anti-instanton pair. Fig. a shows logdet(^) as a function of z (in units of p) for cos 9 = 1, 
Fig. b gives the dependence on cos 6 for z = 1, and Fig. c the behavior of the determinant as a 
function of cos0 for cos 6 = 1/2 and z = 1. 

FIG. 3. Spectrum of the fundamental and adjoint Dirac operator in a random instanton 
ensemble. The spectral density is given in arbitrary units. 

FIG. 4. Quark condensate in an interacting instanton ensemble as a function of quark or gluino 
mass. Fig. a shows the quark condensate for Nf = 1, ... ,4 Dirac fermions in the fundamental 
representations SU(2). Fig. b shows the gluino condensate for N a d = 1,2 Majorana fermions in 
the adjoint representation. 

FIG. 5. Spectrum of the fundamental and adjoint Dirac operator in an unquenched instanton 
ensemble. In the case of the fundamental spectrum the ensemble was created with Nf = 2 fun- 
damental Dirac fermions, while in the case of the adjoint spectrum the ensemble corresponds to 
N a d = 1 adjoint Majorana fermions. 
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FIG. 1. 
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FIG. 2. 
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